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ABSTRACT 


The feasibility of using a miniature heater for a zinc oxide 
thin film oxygen partial pressure sensor was investigated. A silicon 
heater 0.010 inch square by l/2 inch long was used to obtain a t00° C 
surface temperature. The temperature of the heater can be measured by 
the change in resistance of the heater at its operating temperature. 
Experimental data are presented for the re si stance -temperature charac- 
teristics, the temperature distribution, and power requirements of the 
heater. This work has shown that this heater concept is feasible, but 
more effort will be required to deposit the zinc oxide film on the 


heater surface. 
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CHAPTER I 


INTRODUCTION 

2.1 of our past space flights have been flown with a nearly 
100 per cent oxygen atmosphere. This atmosphere was selected for the 
Mercury flights because it allowed the atmospheric supply and control 
system to have the lowest possible weight. The small thrust capabilities 
of the first manned launch vehicles made this a very important decision. 

I. The Need for Oxygen Sensing 

As the length of the space mission is extended, physiological and 
fire hazard considerations become more important than the weight savings 
of a nearly 100 per cent oxygen atmosphere. There are a number of con- 
flicting problems in selecting a cabin atmosphere. The cabin pressure 
should be an absolute minimum to reduce the pressure difference across 
the pressure hull. This would allow the pressure hull to have the mini- 
mum possible weight. For the health of the astronauts in a long-term 
mission, the total pressure and oxygen partial pressure must be within 
certain limits. Also, the fire hazard is increased as the oxygen partial 
pressure is increased. 

The atmosphere selected for the Advanced Apollo Applications and 
Manned Orbiting Laboratory programs consists of a 5 psia total pressure 
with about 65 per cent oxygen. AAP uses a nitrogen diluent and MOL uses 
a helium diluent. 

Current extravehicular activity and lunar excursion space suits 
use a pure oxygen environment at 3 psia. This atmosphere was selected 
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mainly because mobility of current suits is greatly reduced at higher 
pressures. Once this problem is eliminated, all current indications 
point toward the use of a reduced oxygen atmosphere for this application 
also. 

Any atmosphere other than 100 per cent oxygen requires an oxygen 
sensor in order to measure and control the required level of oxygen. 

II. Problems in Sensing Oxygen 

When the requirements of space flight are applied to the possible 
methods of accurately sensing oxygen partial pressure, a very limited 
number of currently known methods remain feasible. The sensing and con- 
trol system must have low weight, small size, and smell power consumption 
plus the ability to operate in a zero gravity condition and withstand the 
vibration and acceleration loads of flight. In addition, the sensor must 
have an electrical signal output that is specific only to oxygen with at 
least a 2 to 5 per cent accuracy and stability over a 100-day period. 

These requirements rule out most of the standard chemical analysis 
techniques that might at first be considered and leave techniques based 
more on the physical and electrical properties of oxygen. 

III. Available Sensing Techniques 

* 

The best commercially available laboratory standard oxygen sensor 
makes use of the paramagnetic properties of oxygen. Not only is oxygen 
strongly paramagnetic, but also, no other gas with paramagnetic or dia- 
magnetic properties is present in sufficient quantities to influence the 
sensor output. This technique would be perfect for spacecraft applica- 
tions except that the magnetic forces involved are so small that a 
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delicate measuring setup is required to detect them. A delicate instru- 
ment cannot pass the vibration and acceleration specifications for launch. 

The polarographic oxygen sensor makes use of an electrochemical 
reaction to produce an output proportional to the cxygen partial pressure. 
While this sensor meets most of the flight requirements, it is subject to 
considerable output drift over long time periods. 

Oxygen has a usable absorption frequency in the ultraviolet range. 
A considerable amount of effort has been expended to design a flight 
sensor using this technique, but no workable device has been produced. 

Another sensor, termed a solid electrolyte sensor, makes use of 
a technique involving the conduction of ionic oxygen through a ceramic 
crystal lattice. A differential voltage is produced that is proportional 
to the ratio of the oxygen partial pressures on the two sides of the 
ceramic. In order to achieve reasonable stability and freedom from 
poisoning, this sensor must be operated at 850° Centigrade. Although a 
commercial sensor using this technique is available, efforts to minia- 
turize the sensor and reduce its electrical power requirements have led 
to many problems achieving hermetic seals at the 85O 0 operating 
temperature. 

A very successful sensor, which is capable of being flight quali- 
fied, has been developed using a magnetic sector mass' spectroscopy tech- 
nique. This sensor measures oxygen, nitrogen, carbon dioxide, and water 
vapor partial pressures. It weighs approximately 5 pounds and requires 
about 5 watts of power. This sensor is very well suited for sensing a 
spacecraft atmosphere, but other applications, such as future lunar 
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excursion and extravehicular activity space suits, will require an oxygen 
sensor that is considerably lighter and less demanding of power. 

The sensor which seems to have the greatest promise in being able 
to successfully meet all requirements is the zinc oxide thin film oxygen 
sensor. 


XV. Zinc Oxide Thin Film Oxygen Partial Pressure Sensor 

Zinc oxide has been used for many years as a pigment for paint, 
as an activator for the vulcanization of rubber, and as a catalyst and 
photocatalyst for chemical reactions (Ref. 1, p. 193)* A considerable 
amount of research has been done on this material in the last 30 years. 

In 1950, Bevan and Anderson (Ref. 2) reported data on the variation of 
electrical conductance of zinc oxide with oxygen pressure at different 
temperatures. Kefeli in 1957 (Ref. 3) proposed using this effect to 
measure oxygen. 

Under Contract NASl-7087, Research Triangle Institute conducted 
an investigation of zinc oxide thin film for sensing oxygen and con- 
structed a laboratory model using this technique (Ref. ^). Their working 
model consisted of a nichrome wire heating element surrounded by a hollow 
cylindrical substrate on * ich a zinc film was deposited and then oxidized 
at 600° C. This whole unit was then enclosed in an insulating housing 
similar to a miniature thermos bottle. This arrangement resulted in a 
sensor weighing about 30 grams, occupying approximately 30 cubic centi- 
meters of volume and requiring 5 to 6 watts of power for operation at 
400° c. 
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In general, the zinc oxide operates as a doped semiconductor. A 
number of different models can each be used to describe the properties 
of zinc oxide under different conditions. If the oxygen is considered 
to diffuse into the bulk material or be adsorbed on the surface, the end 
result is to form new acceptor sites and thus reduce the density of 
conduction electrons in the n-type film. Since the conductivity is pro- 
portional to the electron density, the conductivity is inversely propor- 
tional to the oxygen partial pressure. A much more detailed treatment 
of the theoretical operation of zinc oxide thin films is given by Royal, 
Wortman, and Monteith (Ref. 4 , pp. I7-3O). 

The electrical conductivity and the time required to respond to 
a change in oxygen partial pressure are both temperature dependent. The 
response time decreases with increasing temperature. The temperature 
dependence of the electrical conductivity is shown in Figure 1. The 
minimum in the curve at 400° C is the best operating point for this tech- 
nique for two reasons. First, the time required for the sensor output to 
respond to 50 per cent of its final value for a change from pure oxygen 
to pure nitrogen atmosphere varies from 30 seconds to 6 minutes at this 
temperature, depending on the direction of the change. Secondly, the 
slope of the conductance versus temperature curve at this point is small. 
As can be seen in Figure 1, 60° and 100° temperature changes result in 
only 1 and 5 per cent changes, respectively, in electrical conductance. 
This relatively weak dependence upon temperature at this point minimizes 
the required temperature control accuracy and the effect of a temperature 
gradient along the required heater (Ref. 4, pp. 43-44). 



Conductance (mho) 



Figure 1. - Electrical conductance of zinc oxide vs increasing temperature 
for various partial pressures of oxygen. (Ref. k, p. 5) 
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V. Concept of Miniaturised Heater for 
Thin Film Sensor 

In any design of a high -temperature heater, the required heater 
input power is inversely proportional to the amount of insulation present 
and proportional to the surface area of the heater* Minimizing the 
surface area will minimize the convective and radiant heat losses. 

The convective heat loss is given by 

1 0 - A s h c< T - V (1) 


where 

A g = surface area 

h c * film coefficient 

T * surface temperature in °K 

T 0 = surrounding air temperature in °K 

The radiant heat loss is given by 


= A sM T * T o> 


( 2 ) 


where 



a 


= emissivity 


(Ref. 5* PP» 351*202-211) 


In decreasing the cross-sectional dissensions of the heater in order 
to decrease the surface area, two problems become evident. The heater 
material will have to be selected so that it will not require excessive 
driving voltage or current. High voltages or large currents could only 
be provided by a special converter whose weight and volume would be charged 
to the sensor. Also, small size aggravates the problems of electrical 
connections and heater temperature sensing. 

The objectives of this work were as follows: 

1. Minimize heater power by making dimensions as 
small as practical. 

2. Select a heater material with proper resistivity. 

3* Develop reliable electrical connections to the heater. 

4. Develop a method of measuring the heater temperature. 
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SELECTION OP HEATER MATERIAL 

I. Approximate Range of Resistivity 

The following parameters were chosen as the initial design criteria 
for a heater consisting of a single short length of conductor: (1) volt- 
age less than 100 volts, (2) current between 1 and 100 mill! amperes, 

(3) resistance between 300 and 5*000 ohms, (l*) length of the heater 
between l /4 and 1/2 inch, and (5) heater diameter between 0.001 and 0.050 
of an inch. 

Using these values, the maximum acceptable resistivity is 100 ohm- 
centimeters and the minimum is 0.002 ohm-centimeters at 400 ° C, a range 
falling within the properties of semiconducting materials. Since silicon 
is probably the most readily available and most highly documented semi- 
conducting material, it was selected for the heater. 

II. Silicon Heater 

Since the resistivity of any semiconductor is a function of its 
doping level in the extrinsic region, the resistivity can be adjusted to 
achieve the desired miniature sized heater. 

The effect of temperature on the resistivity of silicon for various 
doping levels is shown in Figure 2 . Silicon has a positive temperature 
coefficient of resistivity in its extrinsic range. In the intrinsic 
region the temperature coefficient of resistivity is negative and the 
logarithm of the resistivity is directly proportional to the reciprocal 
of the temperature. 
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The room temperature resistivity of silicon can be selected so 
that at 4 00° C the silicon bar is operating in its intrinsic region. A 
cross-sectional dimension of 0.010 inch was selected to simplify fabrica- 
tion problems. 

In order to make use of the constant tenperature coefficient of 
log resistance for the intrinsic region to sense the temperature, as 
discussed in Chapter 3V, a room temperature resistivity of approximately 
1 ohm- centimeter was chosen. With a heater length of l/2 inch, and 
assuming the whole bar to be at 400° C, the heater resistance should be 
approximately 2,500 ohms. 

A silicon bar with a room temperature resistivity of 0.8 ohm- 
centimeter and the following dimensions was procured. 

o.et" 0 . 0 /" 
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CHAPTER III 

HEATER ELECTRICAL CONTACTS 
I. Achieving Electrical Contact 

Gold wires O.OOl inch in diameter were first used in an attempt to 
make electrical connections to the bar. They were thermocompression 
bonded to an aluminum pad deposited on the bar ends. Ibis is a standard 
technique used in microelectronic circuitry, and the contacts were 
expected to be capable of carrying the necessary current. However, the 
wires melted just before the thermocompression bond bead at a current of 
11 milliamperes. An attempt to parallel three 0.001-inch gold wires was 
also unsuccessful; all three melted at a current of 18.7 milliamperes. 

Liquid bright gold paste leads fired on the bar proved to be a 
successful method of achieving electrical contacts. The gold paste was 
painted on the ends of the bar, down the ceramic supports, and on to 
conducting pads on the base. 

Liquid bright gold contacts have been known to degrade in a few 
days when operated at elevated teng>eratures (Ref. k, p. 2)» In this 
particular application no degradation has been noticed over periods as 
long as 6 months. The ends of the bar where the contacts are located 
are not at a temperature anywhere near l^OO 0 C. Ibis is the most likely 
reason for the stability. 

Although the electrical contact problem was solved, the bar resis- 
tance was found to be a function of the current for currents small enough 
to cause no measurable temperature rise. 
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II. Ohmic Contact Problems 

When a junction is formed between a metal and a semiconductor, 
either one of two types of contacts are formed. 

If the work function of the metal is greater than the work function 
of the n-type semiconductor, electrons flow from the semiconductor to the 
metal when contact is made. This flow continues until a negative charge 
is built up and the Fermi levels are continuous across the junction. This 
forms a potential barrier at the junction which can be increased or 
decreased by applying an external potential of the proper polarity. The 
magnitude of the net current flow is inversely proportional to the barrier 
potential which is influenced by the polarity of the applied voltage. 

Thus a rectifying junction is formed. 

If the work function of the metal is less than the work function 
of the n-type semiconductor, electrons flow from the metal to the semi- 
conductor. This results in the lack of any potential barrier at the 
junction. The magnitude of the net current flow is independent of the 
polarity of the applied voltage. This is characteristic of an ohmic con- 
tact between the metal and semiconductor. 

In actual practice the height of the potential barrier is more 
dependent on the surface states of the semiconductor than on the metal 
work function. If the surface state density is large enough, contact 
with the metal causes electron flow between the metal and the surface 
states instead of the bulk material. If tip 3 surface states are not 
removed, a rectifying contact will be formed (Ref. 7, pp. 227-235)* 
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When electrical contact Is made to each end of a semiconducting 
material with rectifying contacts, one or the other of the two rectifying 
contacts will he reverse biased for either polarity of applied voltage. 

The current flowing in each direction across a metal semiconductor 
junction is dependent on the thermionic emission of electrons across a 
potential barrier. This current is given by the Richardson -Dushman 
equation 

J = A q T 2 e _E w/ kT (3) 

where 

J ■ current density 
A 0 * 1.2 X 10^ amps/m 2 -°K 
T = temperature 

E w = height of potential barrier 

k = Boltzmann ! s constant (Ref. 7> PP* 232,313) 

When an electric field is applied, the thermionic emission is 
influenced by the space charge of the emitted electrons. This is called 
Schottky effect and is given by 

J * A Q T^ e“GV" c ( c £/ 1+716 ) ^ J/ kT (4) 

where 


c = electron charge 
£,= electric field 
e = permittivity 


(Ref. 7, p. 317) 
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In a rectifying contact, the current from the metal has to over- 
come a potential harrier equal to the difference between the metal and 
the semiconductor work functions. This current is then 

I^g = BAqT 2 e"C E ™' Es ' K ^ 2 3/ kT (5) 


where 

B = area of the contact 
% = metal work function 
E s = semiconductor work function 



The current from the semiconductor has to overcome the barrier 
equal to - E s - cV> where V is the applied voltage. V is positive 
for forward bias and negative for reverse bias. This current is then 

I s _^ m = BA q T 2 / kT (6) 

The net current is then 


1 = ls-*m * lmr»s 

Since one of the contacts will always be reverse biased, it will 
limit the current for the whole bar. For this situation, the forward 
current Ig-^ will be negligible compared to the reverse current. The 


net current will be 
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I = BAjjT 2 e '& n " Es " K ^ ^ ^ kT 

*“ I - InfBAoT 2 ]- F* j^ ~ s) + 

Assume 



kT 


Then 

IB I = K’E 1 / 2 + c (7) 

The characteristics of the contacts for the 0.01-inch silicon bar 
with fired-on liquid bright gold contacts is shown in Figure 3* This is 
a plot of the bar current as a function of the square root of the voltage 
applied to the bar for currents sufficiently small so as not to cause any 
measurable rise in the bar temperature. For these low currents, the 
voltage drop across the bar can be neglected compared to the drop across 
the reverse biased rectifying contact. If it is assumed that £= V/x, 
where x is the width of the space charge region, equation (j) would 
then become 

In I = K’x* 1 /^ 1 / 2 + C (8) 

Although both curves in Figure 3 have straight line segments, the 
slope of the lines increases at 1. 5 to 2.0 volts^/^. This could be 
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accounted for by a decrease in the space charge width as the voltage 
is increased. 

These data are sufficiently close to the characteristics of a 
rectifying contact to verify that the contacts are responsible for this 
effect. 

As the current and therefore the bar temperature increase, the 
rectifying contact will break down to the point where it has a negligible 
effect on the bar resistance. At the 400° operating temperature, the 
resistance characteristics will be due to the silicon bar alone. 



CHAPTER IV 


MEASUREMENT OP HEATER TEMPERATURE 

The temperature of a heater as small as 0.01 inch is very hard 
to measure by conventional means without establishing an undesirably 
large heat drain. 


I. Change in Resistance of Heater to Sense 
Heater Temperature 

As shown in Figure 1, page 6, the logarithm of the resistivity of 

silicon is a linear function of inverse temperature for the intrinsic 

region. Figure 4 shows the resistivity of silicon plotted against 

temperature instead of inverse temperature. The change in the slope of 

o 

the intrinsic curve from about 350 to 450 C is small. Thus, over this 
limited range the logarithm of the resistivity can be assumed to be 
linearly dependent on temperature and the logarithm of the conductivity 
win be directly proportional to temperature. 

Thus a signal proportional to the temperature of the heater in 
the intrinsic region can be obtained by dividing a signal proportional 
to the heater current by a signal proportional to the heater voltage. 

A standard analog computational device is commercially available whose 
output is proportional to the logarithm of the ratio between two input 
signals. 

Once the re si stance -temperature characteristics of the heater are 
known, this latter method would allow the heater temperature to be meas- 
ured completely external to the heater. Although there will be a 
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Figure k,* Resistivi 
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temperature gradient along the heater, the relative insensitivity of 
the zinc oxide film to temperature variations at 400° C would minimize 
this problem. The actual temperature distribution for the heater is 
discussed in Chapter VI. 

II. Practical Experimental Temperature Measurement 

Although the heater temperature will be measured in use by the 
foregoing resistance change method, this method requires a calibration 
technique for determining the bar re si stance -tenperature characteristics. 

One of the first methods considered for measuring the temperature 
was the use of an infrared spectrophotometer. A survey of commercially 
available instruments was made. Although it was possible to focus such 
an instrument down to 0*.Q10 inch, a temperature of 400° C was below its 
accurate usable range. 

An attempt was made to heat the bar in a furnace and measure its 
resistance as a function of the furnace air temperature. Due to the 
rectifying contacts which caused the measured resistance to be a function 
of the measuring current, this method was discarded as being unreliable. 
Also, since in operation the whole bar would not be at the same tempera- 
ture, this method would not give directly applicable results. 

The method finally selected consisted of using a commercially 
available microminiature copper constant an thermocouple Junction. This 
thermocouple consisted of 0.001-inch-diameter wires sheathed in a 
0.030-inch-diameter stainless steel tube. The wires protrude beyond 
the tube and are formed into a Junction of 0.010 inch diameter. The 
thermocouple was rigidly mounted, and the heater bar with its ceramic 
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holder was mounted on the table of a three -axis micropositioner. With 
the aid of a 120 power stereo zoom microscope, the thermocouple was 
placed ©gainst the surface of the heater. TOiis was measuring the sur- 
face temperature which is the temperature of interest for the zinc oxide 
film. 

This is not an ideal way of measuring the temperature, since the 
thermocouple is about the same size as the heater. The positioning of 
the themocouple was critical since any surface irregularities would 
influence its readings. 

Placing the thermocouple on the bar surface would also conduct 
a certain amount of heat away from the bar and, thus, lower the surface 
temperature. This heat conduction appeared to be small enough not to 
influence the data since no change in the heater voltage or current was 
noted when the thermocouple was placed on the surface. 



CHAPTER V 


TEMPERATURE CHARACTERISTICS OF SILICON 
The conductivity of a semiconductor is given by 

o - no(i e + poMjj (9) 

where 

c - conductivity 

« 

n * electron density 
c * electronic charge 
H e * electron mobility 
p * hole density 

= hole mobility (Ref. 7, p. 203) 

I. Ofoeoretical Discussion of Conductivity 
Versus Temperature 

The n-type silicon bar used for the heater had a room temperature 
resistivity of approximately 0.8 ohm-centimeter. This requires a doping 
level of about 7.0 x lO 1 ^ iu^purities/cm^ (Ref. 6 , p. 19 )• 

In the extrinsic region, the electron mobility is an extremely 
complex function of the impurity concentration and temperature for 
impurity concentrations as large as lO 1 ^ irapurities/cm^. Theory cannot 
satisfactorily predict this dependence (Ref. 6, p. 16). The most 
reliable source f data in this region is experimental data such as is 
given in Figure 2, page 10, but actually this consideration is only of 
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limited importance for the operation of the heater since it will not he 
used in this region. 

In the intrinsic region the electron and hole densities are equal 
and the Fermi level lies midway between the valence and conduction bands. 
The electron and hole concentrations are given by 


n 


\ 

N 0 e- (E c- E f>A T 


( 10 ) 




( 11 ) 


where 

N c = effective density of states in conduction band 
N v * effective density of states in valence band 
E c = energy of conduction band edge 
E y * energy of valence band edge 

Ef * Fermi energy (Ref. 7, pp. 199>200) 

Assume 




where 

Eg * width of the forbidden energy gap 


Since 


E f * 


(Eg + Ey) 
2 


E C ■ E f “ E f ■ E v = -£ 


Then 
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Since 



a 


p - F_1 r>V 2M 

|N C e(Pe + p h )J 

in p = -In(N c e(p e + p,)] + 3l (X2) 

The change in the first term with temperature is small compared, 
to the second term, so that a plot of In p versus l/T will be a 
straight line of slope Eg/2 k. 

II. Resistance Versus Temperature 
Experimental Data 

Figure 5 shows the experimentally measured dependence of the 
logarithm of the resistance of the silicon bar heater on temperature. 

The 0.010 -inch- square by 1/2-inch-long bar was mounted in an alumina 
frame with fired-on liquid bright gold electrical contacts. A variable 
constant current power supply was used to supply power to the heater. 

The temperature was measured by a microminiature copper constant an 
thermocouple as described in Chapter IV. Two separate units were tested 
under the same conditions. 

The initial sharp decrease in resistance in Figure 5 at room 
temperature is due to the rectifying contacts on the silicon heater as 
discussed in Chapter III. Once the rectifying contact breaks down, the 
extrinsic curve for both heaters takes on the shape comparable to the 
extrinsic region in Figure k, page 20. 
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The transfer from extrinsic conduction occurs at approximately 
200° C, which is the theoretical value for silicon (Ref. 7, p. 209). 

The intrinsic curve for both heaters is a straight line with a 
negative slope., This is comparable to the similar region of Figure 4, 
page 20. The slopes of the two curves are approximately the same, but 
one heater bar has a higher resistance than the other. This is most 
likely due to physical differences between the two heaters. The silicon 
bars were cut from a 0.0l6-inch-thick wafer and then chemically etched 
down to a nominal 0.010-inch- square cross section. A 20 per cent 
difference in the cross-sectional dimension would be sufficient to cause 
this difference in resistance. When the possible variation in length 
and cross-sectional dimensions are considered, these differences are 
reasonable. 

Equation (12) shows that the slope of a curve of Zn p versus 
l/T in the intrinsic region should be equal to Eg/2k. The dashed line 
in Figure 6 is a representative slope that a sample of silicon with an 
energy gap equal to 1.1 electron volts should exhibit. The two plotted 
curves are the experimental data from Figure page 26, plotted against 
reciprocal temperature. The theoretical slope is 


31ope theoretical 


i — = 6.58 x 10 3 

(2) (8. 617 x 10-5) 


The slope of the experimental data is 


sl0 P e exp e rimental 


(Zn 10^ - Zn 1.2 X 10 3 ) = ^ x 1Q 3 

(2.18 X lO' 5 - 1.5 x 10"5) 



In Resistance 
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This difference in slope is due to the fact that in the experi- 
mental case the heater bar is not all at the same temperature. The 
experimented curve is for the average resistance for the whole heater 
plotted against the maximum temperature. 

The mathematical model developed in Chapter VI was used to 
calculate the temperature distribution along the heater for three 
different maximum temperatures. Using the theoretical variation of 
resistivity with temperature for the intrinsic region and two straight 
line approximations of the data in Figure 2, page 10, for the extrinsic 
region, the average resistivity of the bar was calculated. The slope 
of the mathematical model data plotted in Figure 6, page 28, is slightly 
smaller than that of the experimental data. The mathematical model 
slope is 

Slope model “ 5>0 x 103 

This is only 3.6 per cent smaller than the experimental slope. This is 
sufficiently close to verify that the slope of the experimental data is 
due to the temperature distribution along the heater. 


CHAPTER VI 


TEMPERATURE DISTRIBUTION ON THE HEATER 


As discussed in the last chapter, the temperature distribution 
on the heater bar has a significant effect on its operation. The temp- 
erature distribution will be discussed on the basis of a mathematic 
model and experimental data. 

The following assumptions have been made for this model: 

1. No temperature variation through heater 
cross section. 

2. Heater is symmetrical about its center. 

3* Maximum temperature exists at the center. 

4. To avoid trial and error solution, the 
boundary conditions will be taken from the 
experimental data. 

5. Heat flow into a volume element is positive. 

The heat balance equation for any volume element in the heater 
consists of 


q x + q 2 + qg + q r + <l c = 0 


(13) 


where 


*2 

’g 


= rate of conductive heat flow into the volume element 
in Btu/hr 

= rate of conductive heat flow out of the volume element 
in Btu/hr 

= rate of heat generation inside the volume element in Btu/hr 
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q r * rate of radiant heat loss from the surface in Btu/hr 

q c = rate of free convective heat loss from the surface 
in Btu/hr 


,l=-kA^ 

1 dx x 


%2 “ +kA 


dT 


= kA — I + kA 


^ *x+Ax ^ *x 


d 2 T 

dx? 


qg - - 3 -411^ 


2„ dx 


Prom equation (2), 


q_ = 4Dh_(T - T 0 )dx 


From equation (1), 


% = -4Dh c (T - T 0 )dx 

Combining equations (13), (14), (15 ) , (16), and (17) gives 


kA 


d 2 T 

dx 2 


3-4lI 2 P 
+ A 


- 4D(h r + h 0 )(T - T 0 ) = 0 


(14) 

(15) 


(16) 


(17) 


(18) 


The quantities k, p, h r , and h c are all temperature dependent so 
equation (l8) is a nonlinear differential equation. The most practical 
solution for this equation is a numerical solution with the following 
approximations for the temperature dependent factors: 



For T > J00° C 
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p _ . 2.49 X 10- 6 e 6 *3 8xlo5 /( T+2 75) 

ohm-ft 

For 200° C < T < 300° C 


p ohm-ft = 0*383 e 


60x10"3t 


For T < 200° C 


p s. o. 07*54 ^ 

p ohm-ft ,,0T54 e 

The approximation for T > 500° C is from the intrinsic curve of 
Figure 2, page 10, and the other two are from the extrinsic curve in 
Figure 4, page 20. 


k Btu/hr-ft-°F 


3.38 X 10 5 t" 0 , 5T 

(Ref. 8 , p. 99) 


From equation (2) 



oF t 


where 



(Ref. 9, p. 70) 


(Ref. 5, p. 314) 
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The computer program and the data readout for the numerical 
solution for the mathematical model are given in the appendix. Figure 7 
is a plot of the mathematical model and experimental temperature distri- 
bution. The experimental distribution agrees reasonably well with the 
model distribution for the middle region of the bar. The end region 
temperatures are of little importance for the operation of the sensor 
since in order to keep the temperature variation along the zinc oxide 
film within approximately 50° C> the film would be deposited only within 
0.050 inch of the center. Within this region the model is sufficiently 
accurate to predict the temperature distribution. 






CHAPTER VII 


POWER REQUIREMENTS OF HEATER 


Because of its small size, the heater car. he operated in a com- 
pletely uninsulated condition. Since silicon has a negative temperature 
coefficient of resistance in the intrinsic state, the heater must he 
operated from a constant current power supply. If operated from a 
constant voltage supply, a completely unstable condition will result 
at temperatures about 200° C and the temperature will increase quickly 
to the melting point for silicon. 

Figure 8 shows the current and corresponding voltage required 
as a function of the maximum heater temperature for a typical heater. 

A 50-milliamp constant current power supply capable of putting out 

100 volts would be required to power the heater with no forced convection 

and no insulation. 

Figure 9 shows the power required to heat the heater in the 
uninsulated condition as a function of maximum temperature. The heater 
required approximately 2.5 watts of power to reach a maximum temperature 
of k 00° C. 

In a final configuration this power requirement should be able 
to be reduced through the use of a reflecting radiation shield around 
the heater. 
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Current 



Figure 8. - Current -voltage requirements of the heater 





CHAPTER VIII 


CONCLUSIONS 


This work has shown that the use of a small-diameter heater for 
a zinc oxide oxygen sensor is feasible and will result in a considerable 
reduction in size and power requirements. 

The temperature of the heater can be sensed by measuring its 
average resistance. A significant variation in the resistance charac- 
teristics for different elements has been noted, but the slope of the 
log resistance versus temperature curve for the operating temperature 
region is repeatable between different heaters. Also, the maximum 
voltage occurs at approximately the start of the intrinsic region at 
200° C. Urns, if the current is increased from the point of maximum 
voltage until the logarithm of the resistance has decreased by the slope 
of the experimental data divided by the temperature change, an approxi- 
mate operating point can be determined. 

Thus, from the data of Figure 5 , page 26, 

In R » K • 8.12 X 10"5 T (°C) 

A In R = A(-8.12 X 10 ' 5 T) 

R 2 = Bj. e - 8 - 12 XlO-5(T 2 -Ti) 

so that 

R 4C0° C = R 200° C 
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Using the two values of Rgoo 0 C from F:t S ure page 26, of 5.7 K 
and 4.1 K, R^qqo c calculated from equation (19) is 1.12 K and 0.8l K, 
respectively. These values would give a temperature error of about 
10° C. 

This effort has shown mainly that this type of heater is feasible, 
but considerably more work needs to be accomplished before a prototype 
sensor can be built. Eie techniques of electrically isolating the zinc 
oxide film from the heater and depositing the film on the surface still 


have to be worked out. 



REFERENCES 


1. Heiland, G. ; Mollwo, E. ; Stockman, F. ; ’'Electronic Processes in 

Zinc Oxide.” Solid State Physics , vol. 8, Academic Press, 

New York, 1959* 

2. Beven, D. J. M.$ Anderson, J. S. : Discussions Faraday Society . 

8:238, 1950, cited by Ref. 1, p. 2^9* ~ 

3. Kefeli, A.: Diplomarbeit . University of Erlanger, 1957> cited 

by Ref. 1, p. 277. 

4. Royal, T. M. ; Wortman, J. J. \ Monteith, L. K. : "An Investigation 

of Thin Film Oxygen Partial Pressure Sensors. ” NASA CR 1182, 

1968. 

5. Kreith, F. : Principles of Heat Transfer . International Textbook 

Company, Scranton, 195^7 ’ 


6. Integrated Silicon Device Technology Volume I - Resistance . 

USAF ASD-TDR-63-3I6, 19^57 

7. Azaroff, L. V. ; Brophy, J. J. : Electro nic Processes in Materials . 

McGraw-Hill Book Company, New York, 1955. 

8. Integrated Silicon Device Technology Volume V . USAF ASD-TDR-63-516* 

1965 . 

9. International Journal of Heat and Mass Transfer Volume 

Pergamon Press, New York, 1952V 


ko 


41 


JOB NO. 1302 ASSIGNED TO JOR, 1, 100, 4 0000. A1065 2, CARL R PEARSON T 

START OF OUTPUT FOR JOB NO. 1302. 

TER1302. 

JOB, 1,100, 40000. A 1 0 0 5 2, CARL R PPAPSO^' TBC125,120? 

PROGRAM TEMP( | w PUT, OUTPUT) 

000003 DIMENSION Y (200),T(200 ), X (200) 

000003 TAMD-75. 

000004 T(l)-743. 

000006 0X-4.15P-O4 

000007 AR“B , 9 5F-07 

000011 0*8 , 34E-04 

000012 Cl 043 

000014 D0300 I *1, 3 

000016 RUQT*0. 

000017 Y(l)-". 

000017 X(l)-0. 

000020 TC«(T(l)-32. )*5./9. 

0000 24 P1IHT2,TAMtt,T(l>,TC,)C<l> 

000040 2 FORMAT(//13X, 12MTEMO. DEG , F,6X, 12MTFMP. RFC, C,5X,15H»»0SIT|0« I M C 

lMES,5X,OHTAMBIEMT«Ein.3/14’<,F10.3,8X,Ein.3,OX,Ein.3) 

000040 Hl-41 

000041 DO 10 ON *2, Ml 

000043 X(N)«X(N-l)+nx*12. 

000046 A*T(N-1) 

000047 B*A-TAMR 

000051 HG-.72*.172*C((A+460. )/100. ) **4-019. 25)/^ 

000060 HC“.27*(n/D)**,25 

000065 TK=3.38 p 03/A**,57 

000071 IF(A-572. )10,10,30 

000074 10 IF(A-302. >11,11,12 

000077 11 RMO-.0764*EXP(5.58E-03*5.*(A-32. )/0. ) 

000107 GO TO 40 

000107 12 RHO*, 3 83 *EXP( -2.60E-03*5,* (A- 32. )/9. ) 

000117 GO TO 40 

0 00117 30 RHO*2.49E-0n*EXP(6.38F03/(<A-3?. )*5,/0. + 2 73. )) 

000130 40 RMOT*RMOT+RHO 

000132 SLY*(4.*D*(HR+MC)*R-3.41*CI **2*R'!D/^R)/ (TK*AR) 

000146 Y ( N ) “ Y ( N- 1 ) +$ LY *D X 

000152 T(N)»T(N-1)+Y(M)*DX 

000155 IF(T(N).LE.TAMR)f'l*M 

000161 TC* ( T ( N )- 32 . ) * 5 , / 9 . 

000165 100 PRI NT1,T(N ),TC, X (N) , R'lOT 

000204 1 FORMAT(14X,Pin.3,0X, c l0.3,OX,E10.3, f )X,F10.3) 

000204 200 T(l)*T(l)-50. 

000206 300 C I »C I - , 008 

000212 STOP 



000214 

END 







PROGRAM LENGTH INCLUDING 

I/O 

BUFFERS 





005551 

FUNCTION 

ASSIGNMENTS 







STATEMENT 

ASSIGNMENTS 







1 

000300 2 

- 

000236 10 

- 

0*0074 11 

- 

0*0*77 

12 

000107 30 

- 

000117 4* 

- 

06*130 2*0 

- 

09*704 

BLOCK MAI 

ES AND LENGTHS 







V ARIARLE 

ASSIGNMENTS 







A 

001476 AR 

- 

001460 R 

- 

001477 Cl 

- 

0*147* 

0 

001467 n,< 

- 

001465 HC 

- 

0*1501 t'R 

- 

0*15*0 

| 

001471 N 

- 

001475 NT 

- 

001474 RMO 

- 

001508 

RHOT 

001472 SLY 

- 

001504 T 

- 

00061' 6 TAMB 

- 

*01464 

TC 

001473 TK 

- 

001502 X 

- 

**1154 Y 

- 

000334 

START OF 

CONSTANTS 







000216 
START OF 

TEMPORARIES 







000310 
START OF 
000325 

INDIRPCTS 







UNUSED C0 1PILER SPACE 







002100 









TEMP. DEG. F 


TEMP. DEO , C 

POSITION |MCH6S 

TA’^MENT* 7.5006+0] 


7.430E+02 


3,950F+*2 


0. 




7.420E+02 


3.0446+02 


4.092E-03 


3.6*06-0? 


7.399E+02 


3.933E+02 


9.0846-03 


7.0296-02 


7.368E+02 


3.9166+02 


1.4006-02 


1.0026-01 


7 . 3 2 S r ; Q 2 


3 . 8S3E+02 


1.0076-02 


1.4206-01 


7, 273E+02 


3.8636+02 


2.4966-02 


1.3106-01 


7.2O8E+02 


7 , R27E + 02 


2.005E-02 


2.7076-01 


7.132E+02 


'-,7346+02 


3.4046-02 


2.02*6-01 


7.042E+02 


■5,?54E + 02 


3.094E-02 


3.0736-01 


G.933E+02 


% , 176E+02 


4.4936-02 


3, 554E-01 


6.813E+02 


1 "106+02 


4.0026-03 


4. *3*6-01 


6.682E+02 


S. 5346+02 


5, 40 IE- 02 


4.E64E-01 


6. 527E+02 


3.443E+*2 


5.0006-02 


5. 3246-01 


6.340E+02 


3.3506+o ? 


5.40*6-02 


5. *? 4 6-01 


6.147E+02 


3.237E+02 


6.0806-02 


5 . 9 3 3 E - 0 1 


5.914E+02 


3. 108 F. + 02 


7.4886-02 


0. *736-01 


5.643E+02 


2.057E+0? 


7.037E-02 


9.4676-01 


S.324E+02 


2.7806+02 


8.4866.-02 


1.124E+00 


4.955E+02 


2.576E+02 


0.086E-*? 


1 .31*6 + 0* 


4, 540E+02 


2.3446+02 


0.485E-02 


1.50*6+00 


4.074E+02 


2. 0856+02 


0.084E-O2 


1.7146+00 


3. 558E+02 


1. 790E+O2 


1,04 66-01 


1.037E+0* 


3.000E+Q2 


1.480E+02 


1.0986-01 


2,1436+0* 


2.409E+02 


1, 161E+02 


1, 14*6-01 


3.316F+00 


1.794E+02 


3.189E+01 


1.1906-01 


2.460E+0* 



1,1625+02 

5.1985+01 


T5MP. nFO. F 
6.9306+02 
G. 9205+02 
6.399E+02 
6.8685+02 
6.8275+02 
6, 7745+02 
6 . 7095+Q2 
6.6325+02 
6.547.5+02 
6. 4385+02 
6. 3185+02 
6 . 1805+02 
6.0235+02 
5.842P+02 
5.6345+02 
5. 3945+02 
5, 1225+02 
4. 8X65*02 
4.4786+02 
4 . 106E+02 
3.7005+02 
3.2645+02 
2.8026+02 
2. 3215+02 
1. 8245+02 
1.316^+02 
7. 9905+01 
2.7725+01 


TEMP. 05P,. F 
6 4305+0? 
6, 4215*02 
6. 4015+02 
6 . 373E+02 
6,3345+02 
6. 2855+02 
6.2255+02 
6.1535+02 
G.Ofi9F+02 
5.9725+02 
5.0605+02 
5.7315+02 
5.5845+02 
5.4185+02 
5.2315+02 


4.079F+01 

1.2485-01 

?. 5796+00 

1.1105+01 

1. 2985-01 

2.6775+00 

IMP, OFH. 0 

*OSITIO»« 1 \'CH*S 

T/kMRIPWT* 7. 

3,6725+0? 

0 . 


3.5675+02 

4 , 9 9?5-03 

6.2976-02 

3.655E+Q2 

9.9845-03 

1.0646-01 

3.6385+02 

1.4985-0? 

1.6085-01 

3.6155+02 

1.9975-02 

2.1675-01 

3.5855+02 

2,4965-0? 

2.747*-01 

3.5505+02 

2.9955-0? 

3.3545-01 

3.5075+02 

3.4945-02 

3.9985-01 

3.457E+02 

3.9945-02 

4,6885-01 

3.3995+02 

4.4935-02 

5.4376-01 

3.3325+02 

4.9925-02 

6.2636-01 

3.2565+02 

5. 4 91*- 02 

7.1905-01 

3.1685+02 

5.99O*-02 

8,2495-01 

3. 0 685+02 

6.4905-02 

9.4915-01 

2.9525+02 

6.9895-02 

1,9995+00 

2. 8195 + 02 

7.4885-02 

1.2765+00 

2.668F+0? 

7.987F-02 

1.4616+00 

2.4985+02 

0.4865-02 

1.65? 6+00 

2.3105+0? 

8.0856-02 

1.8526+00 

2.1035+02 

9.4855-0? 

2,0625+00 

1. 878*+0? 

9.9845-0? 

2.2845+00 

1.6355+02 

1.04S*-01 

2.4995+00 

1.3795+0? 

l.nopp-oi 

?. 6876+00 

1.1115+0? 

1.14?*-01 

2. 8496+00 

8.3555+01 

1. 190*- 01 

2.9896+90 

5.5315+01 

l.?48*-01 

3.1195+00 

2. 6615+ox 

1, 2085-01 

3.712**00. 

'2.3JS5+00 

1.3485-01 

3.3095+00 

:mp, nrn, c 

Pns«T|n^' |momps 

T^M9 | 5)'T = 7 . 

3.3945+02 

0 . 


3,3895+02 

4,9925-03 

8.3255-02 

3.3795+02 

9.9845-03 

1. 0726-01 

3.3635+02 

1.49',5-n? 

2.5285-01 

3.3415+0? 

1.9975-02 

3.4075-01 

3,3145+02 

2.4965-02 

4,3195-01 

3. '>805+02 

2, n°5*-o? 

5.8766-01 

3,2415+0? 

3.4945-0? 

6. '>915-01 

3. 19 4 *+02 

3 . 9945-02 

7.3795-01 

3.1405+0? 

4 . 4955-02 

8.5645-01 

3.0785+02 

4.992*-0') 

9.8726-01 

3.0065+0? 

5.491*-02 

1.1345+00 

2.9255+02 

5 , 9905-02 

1, 303F+QO 

2. 8325+02 

6,4906-02 

i.48?p+nn 

2.7295+02 

6.989*-02 

1. 6655+00 
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5.02 5 c + 02 

2.B14E+0? 

7.4&8E-0? 

1. 8536+1)0 

4 . 7 9 9 F. ♦ 0 2 

?.4f!8* + 0? 

7,np ( 7c.n<> 

P.^'lRC + nn 

4. 552E+02 

2 , 35 1E+92 

?, , 4 o .« C - rt f) 

2. +HPF+nn 

4.235 c +92 

2 . 2036+02 


2. *t p r f + n ° 

3.‘M7F.+92 

2.0436+02 

9 , 4 °* 5 F " ° 2 

2,672.6 + 09 

3 . 6 8 •) r + 0 2 

1.8726+92 

r ', 0 . 0 4 ^ " 7 

7, «?97e + f)0 

3.3626+02 

1.R90F+02 

1.94SF-91 

3.111F+0O 

3.01 JF+02 

1, 4006+02 

l.ORP.F-ni 

3.306^+99 

2.SP2F+02 

1.301E+0? 

1.14TF-01 

3. 4796+09 

2 ,2946+02 

1. 097F + 92 

1.198 c -01 

3,4356+09 

1.917E+92 

8. 4 7 3 6.+ 01 

1 , 2 4 f - n 1 

3. 77'-F+no 

1,532**92 

4.73.5E+01 

1. 29°F» r, l 

3. 

1.1426+0? 

4.5R5F+01 

1.34 "E-oi 

4 , r 1 7 6 + 0 9 

7.466E+ n l 

2.370E+C1 

1. 39°F- , 11 

4,l r '5 r +9 n 


Ob/13/fio lrc cm scop* 3.0 66ooc-i3r< 

10.26. 28. TER1302. 

io.26. 28..j04,i,inn,4nnon. mop? 

10. 20. 23. CARL R P 6 ARSON T9C125,120? 

10. 2U. 29. NOMAP. 

10. 215, 29,ft'J M (S) 

10.26.33. LOO. 

10 , 2t"> . 35 IF'iORY Oil 1 ! OH r.M 
10. 26,36. STOP 

10. 26. 36. CPU 0.502452 SEC . 

10.26. 36. PPU 9.515008 SFC. 


10.2G.36.OATE 06/13/69 


04/14/69 


ENO OF OUTPUT FOR JOR NO. 1302. 
JOR NO. 1302 IS COMPLETE. 



